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Abstract 

Background: Cobalt-ferrite nanoparticles (Co-Fe NPs) are attractive for nanotechnology-based therapies. Thus, 
exploring their effect on viability of seven different cell lines representing different organs of the human body is 
highly important. 

Methods: The toxicological effects of Co-Fe NPs were studied by in-vitro exposure of A549 and NCIH441 cell-lines 
(lung), precision-cut lung slices from rat, HepG2 cell-line (liver), MDCK cell-line (kidney), Caco-2 TC7 cell-line (intestine), 
TK6 (lymphoblasts) and primary mouse dendritic-cells. Toxicity was examined following exposure to Co-Fe NPs in the 
concentration range of 0.05 -1.2 mM for 24 and 72 h, using Alamar blue, MTTand neutral red assays. Changes in 
oxidative stress were determined by a dichlorodihydrofluorescein diacetate based assay. Data analysis and predictive 
modeling of the obtained data sets were executed by employing methods of Knowledge Discovery from Data with 
emphasis on a decision tree model (J48). 

Results: Different dose-response curves of cell viability were obtained for each of the seven cell lines upon exposure to 
Co-Fe NPs. Increase of oxidative stress was induced by Co-Fe NPs and found to be dependent on the cell type. A high 
linear correlation (R 2 =0.97) was found between the toxicity of Co-Fe NPs and the extent of ROS generation following their 
exposure to Co-Fe NPs. The algorithm we applied to model the observed toxicity belongs to a type of supervised 
classifier. The decision tree model yielded the following order with decrease of the ranking parameter: NP concentrations 
(as the most influencing parameter), cell type (possessing the following hierarchy of cell sensitivity towards viability 
decrease: TK6 > Lung slices > NCIH441 > Caco-2 = MDCK > A549 > HepG2 = Dendritic) and time of exposure, where the 
highest-ranking parameter (NP concentration) provides the highest information gain with respect to toxicity. The validity 
of the chosen decision tree model J48 was established by yielding a higher accuracy than that of the well-known "naive 
bayes" classifier. 

Conclusions: The observed correlation between the oxidative stress, caused by the presence of the Co-Fe NPs, with the 
hierarchy of sensitivity of the different cell types towards toxicity, suggests that oxidative stress is one possible mechanism 
for the toxicity of Co-Fe NPs. 
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Background 

Magnetic nanoparticles (NPs) possess unique proper- 
ties which can be applied in nanomedicine: they ad- 
dress targets such as cellular therapy, tissue repair, 
nanobiosensors, drug delivery, magnetic resonance im- 
aging and magnetic fluid hyperthermia. All these appli- 
cations require high magnetization values of NPs and size 
of less than 100 nm with uniform physical and chemical 
properties [1]. Additionally, it is essential to understand the 
biological fate and potential toxicity of magnetic NPs for 
their successful application in nanomedicine [2]. Over the 
years, Iron oxide, especially magnetite (Fe 3 0 4 ), was the 
most investigated magnetic NP. In the last decade it be- 
came easier to synthesize new and more effective types of 
magnetic NPs [3]. Cobalt- ferrite (Co-Fe) NPs which belong 
to the crystal family of spinel ferrites (MFe 2 0 4 ), posses lar- 
ger magnetic anisotropy than other ferrites (e.g. magnetite) 
making them more attractive for nanotechnology based 
therapies [3,4]. 

The toxicity of Co-Fe NPs was explored in a number 
of studies. A significant decrease in cytokinesis-blocked 
proliferation index and increase in the frequency of 
micronucleated binucleated lymphocytes were shown 
when employing 5.6 nm Co-Fe NPs [5]. However, coat- 
ing the surface of these NPs led to a 4-fold reduced level 
of toxicity [5]. Similarly, toxicological study of silica-coated 
Co-Fe possessing a silica shell of 50 nm thickness, revealed 
that although the particles were found in the mice's brain, 
no significant changes in the hematological and clinical bio- 
chemistry tests were found [6]. Another study investigated 
the embryotoxicity of Co-Fe NPs (17 ± 3 nm) through an 
embryonic stem-cell test which show differentiation into 
cardiomyocytes. The obtained ID 50 for the inhibition of dif- 
ferentiation classified Co-Fe NPs coated with gold and si- 
lanes as non-embryotoxic. However, Co-Fe NPs coated 
only with silanes were found to be weakly embryotoxic, but 
less embryotoxic than the cobalt ferrite salt (CoFe 2 0 4 ) [7]. 
Obviously, further toxicological work should proceed with 
the aim to achieve a larger toxicological data-base to enable 
the prediction of toxicology by in-silico approaches [8-10]. 

Predictive toxicology is based on the development of al- 
gorithms that are capable of predicting toxic effects (the 
output) from chemical and biological information (the in- 
put) [11,12]. Knowledge Discovery from Data (KDD) is the 
process of identifying valid, novel, useful and understand- 
able patterns from large or complex datasets. KDD can 
also be applied to small data sets where new insights can 
be inferred, as is the case of this study. The core of the 
KDD process is Data Mining (DM), involving the inferring 
of algorithms that explore the data, develop models and 
discover significant patterns [13]. 

This manuscript explores the toxicological effects of 
Co-Fe (CoFe 2 0 4 ) NPs on viability of cells, representing 
the different organs of the human body, for expanding 



the toxicological knowledge in the future use of bio- 
medical applications. Dose-response curves were car- 
ried out in the concentration range of 0.05 -1.2 mM 
employing MTT, NR and Alamar blue as endpoint as- 
says following exposures for 24 and 72 h. The cell via- 
bility experiments were complemented by determining 
NP-induced changes of oxidative stress in five of the 
cell lines. Finally, we applied KDD and DM to data, 
gathered in the experimental studies, towards the for- 
mation of a predictive model of NP toxicity. The pre- 
dictive modeling was applied by carrying out training 
and validation through an iterative process when apply- 
ing the KDD approach. The model predicts the relative 
hierarchy of the variables studied in the viability tests 
consisting of concentration, cell type and duration of ex- 
posure. A similar approach has been recently reported in a 
study of the toxicity of Co-NPs and Co-ion [14]. The nov- 
elty of the presented model is in its multi-dimensional 
perspective that cannot be achieved by traditional visual 
examination of two or even three dimensional plots, as 
done in most previous toxicological studies. 

Results 

The toxicological effects of Co-Fe NPs were examined 
using seven different cell lines and precision-cut lung slices. 
Since penetration of NPs into the human body proceeds 
principally through inhalation or orally, whereas penetra- 
tion through healthy skin is negligible [15], we have chosen 
cell lines representing lung (A549 and NCIH441 cell lines) 
and intestine (Caco-2/TC7 cell line) as the primary sites of 
interaction. In order to bridge the gap between lung cells 
and the lung organ, we also examined rat precision-cut 
lung slices. Liver (HepG2 cell line), kidney (MDCK cell line) 
and the immunological system (primary murine dendritic 
cells and a human B-lymphocyte cell line - TK6) have been 
selected as the secondary major sites of interaction follow- 
ing the penetration of NPs into the blood circulation. 
Dose-response curves of 0.1-1.2 mM Co-Fe NPs (or 23.5 - 
282 ug/ml) were examined employing MTT, neutral red 
(NR) and Alamar blue as endpoint assays. In addition, ROS 
levels were explored in five of the cell lines. Data analysis 
and modeling of the obtained data sets, for the toxicological 
dose - response curves for Co-Fe NPs were carried out 
using the approach of Knowledge Discovery from Data 
(KDD). 

Cytotoxic effect of Co-Fe NPs on cells 
Cytotoxic effect of Co-Fe NPs on Caco-2 cells 

Caco-2 TC7 cells, incubated with Co-Fe NPs (0.1-1 mM) 
for 24 and 72 h, were analyzed employing three different 
tests for cell viability. These tests, consisting of MTT, NR 
and Alamar blue assays (Figures 1, 2 and 3, respectively), 
show that after 24 h Co-Fe NPs do not significantly influ- 
ence cell viability of Caco-2 cells, however 72 h incubation 
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Figure 1 Effect of Co-Fe NPs on viability assayed by MTT. Cytotoxicity was examined in different cultured cells: (A) Caco-2, (B) HepG2, (C) 
MDCK, (D) A459, (E) NCIH441, (G) TK6 and in lung slices (F). Cells were treated with 0 - 1.2 mM Co-Fe NPs for 24 and 72 h. Cell viability was 
assessed by IvuT assays and results are presented as ratio of control (untreated group). Mean ± SD of three independent experiments performed 
in triplicate. * Significant differences of each point from the control at P<0.05. 



with Co-Fe NPs causes a significant decrease in viability, 
especially at concentration of 1 mM, possessing an average 
EC 50 of 0.86 ±0.02 mM. This decrease is shown by the NR 
and Alamar blue assays and less by the MTT assay. When 
comparing the dose-response curves of the three assays, 
it turns out that the most sensitive one is Alamer blue. 
This assay (Figure 3) shows that Co-Fe NPs at a concen- 
tration of 0.5 mM induce -25% decrease in viability while 



the highest concentration of 1 mM demonstrates viability 
decrease of 50%. 

Cytotoxic effect of Co-Fe NPs on HepG2 cells 

HepG2 cells, incubated with Co-Fe NPs (0.1-1.2 mM) for 
24 and 72 h, were analyzed employing two different 
methods for viability; MTT and NR assays (Figures 1 and 2 
respectively). It is shown that incubation with Co-Fe NPs 



Horev-Azaria et a I. Particle and Fibre Toxicology 2013, 10:32 
http://www.particleandfibretoxicology.eom/content/1 0/1 /32 



Page 4 of 1 7 



O 

o 



o 
o 



Caco-2 



1.2 
1 

0.8 
0.6 
0.4 
0.2 
0 



» 24h 
— B--72h 



0.2 0.4 0.6 0.8 
Co-Ferrite NPs (mM) 



B 

1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 
0 



HepG2 




0.2 0.4 0.6 0.8 1 1.2 
Co-Ferrite NPs (mM) 



MDCK 



1.2 
1 

0.8 
0.6 
0.4 
0.2 
0 



■24h 
- ■- - 72h 



0.2 0.4 0.6 0.8 1 
Co-Ferrite NPs (mM) 



1.2 



Figure 2 Effect of Co-Fe NPs on viability assayed by NR. Cytotoxicity was examined in different cultured cells: (A) Caco-2, HepG2 (B) and (C) 
MDCK. Cells were treated with 0-1.2 mM Co-Fe NPs for 24 and 72 h. Cell Viability was assessed by NR assays and results are presented as ratio 
of control (untreated group). Mean ± SD of three independent experiments performed in triplicate. * Significant differences of each point from 
the control at P<0.05. 



for 24 h does not cause any significant change in cell viabil- 
ity. Incubation with Co-Fe NPs for 72 h reduces cell viabil- 
ity especially at the higher concentration range (1.2 mM). 
The viability decrease is 15% according to the NR assay 
and 50% using the MTT assay. 

Cytotoxic effect of Co-Fe NPs on MDCK cells 

MDCK cells were incubated with Co-Fe NPs (0.1-1.2 
mM) for 24 and 72 h, after which the cells were analyzed 
by MTT and NR assays (Figures 1 and 2, respectively). 
The incubation with Co-Fe NPs for 24 h did not result in 
any significant change in cell viability. However, incuba- 
tion with Co-Fe NPs for 72 h reduced cell viability signifi- 
cantly (50% reduction) at 0.6 mM and 0.8 mM as shown by 
MTT and NR assays, respectively. Cell viability reaches at- 
tenuation of 70-80% reduction at concentration of 1.2 mM 
Co-Fe NPs. The average EC 50 is 0.74 ± 0.19 mM for 72h 
of incubation. 

Cytotoxic effect of Co-Fe NPs on A549 and NCIH441 cells 

A459 and NCIH441 cells were incubated with Co-Fe NPs 
(0.05-1 mM) for 24 and 72 h, and analyzed by the MTT 
assay (Figure 1). The highest concentration of Co-Fe NPs 
(1 mM) causes only 10% and 30% decrease in A549 cell via- 
bility after 24 and 72 h respectively. However, whereas 24 h 



incubation with NCIH441 cells did not cause any signifi- 
cant toxic effect 72 h incubation with Co-Fe NPs caused 
40% and 60% viability decrease at concentrations 0.4 mM 
and 1 mM respectively, possessing EC 50 of 0.51 mM for 
72h of incubation. Thus, the A549 cell line emerges as a 
less sensitive cell line than the other cell lines and especially 
in comparison to the other lung cell line NCIH441 which is 
highly sensitive. 

Cytotoxic effect of Co-Fe NPs on TK6 cells 

TK6 cells were incubated with Co-Fe NPs (0.05-1 mM) 
for 24 and 72 h, and were analyzed by MTT and Alamar 
blue assays (Figures 1 and 3 respectively). This cell line of 
B-lymphocytes seems to have the highest sensitivity for 
Co-Fe NPs, where already at concentration of 0.25 mM 
40% of viability decrease has been found after 24 h incuba- 
tion; after 72 h incubation 0.25 mM cause a viability de- 
crease of 90%. The average EC 50 is 0.37 ± 0.07 mM and 
0.15 ± 0.03 mM for 24 and 72 h of incubation, respectively. 

Cytotoxic effect of Co-Fe NPs on Dendritic cells 

Primary Dendritic cells were incubated with Co-Fe NPs 
(0.05-1 mM) for 24 and 72 h, and analyzed using the 
Alamar blue assay (Figure 3). These cells which also be- 
long to the immune system, as TK6, were found to be 
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Figure 3 Effect of Co-Fe NPs on viability assayed by Alamar blue. Cytotoxicity was examined in different cultured cells: (A) Caco-2, (B) 
primary dendritic Cells and (C) TK6 were treated with 0 - 1 mM Co-Fe NPs for 24 and 72 h on all different cell lines. Cell Viability was assessed by 
Alamar blue assays and results are presented as ratio of control (untreated group). Mean ± SD of three independent experiments performed in 
triplicate. * Significant differences of each point from the control at P<0.05. 



non-sensitive to Co-Fe NPs. The highest concentration 
of 1 mM does not cause any decrease in cell viability, ei- 
ther after 24 h incubationor after 72 h incubation. 

Cytotoxic effect of Co-Fe NPs on lung slices 

Precision-cut lung slices from rat were incubated with 
Co-Fe NPs (0.05-1 mM) for 24 h and were analyzed by 
the WST-1 assay (Figure 1). The viability decrease of 
about 20% already manifests at 0.4 mM Co-Fe NPs. 
These lung slices' sensitivity towards the NPs is demon- 
strated by the lethal effect of an almost 100% decrease in 
cell viability at NP concentration of 0.8 mM. The EC 50 
for precision-cut lung slices is 0.54 mM. 

Leaching of cobalt ions from Co-Fe NPs 

In order to be able to discriminate between the direct 
toxicological effect of Co-Fe NPs and their indirect effect 
due the release of Co-ions arising from the dissolution 
of the NPs, we determined the extent of release of Co- 
ions from Co-Fe NPs. The exact percentage of leaching 
was determined for 1.2 mM Co-Fe NPs (the highest con- 
centration that was used in all toxicological assays). Co- 
Fe NPs were incubated in DMEM with or without 10% 
serum for several minutes and for 72 h. The total 
amount of Co-ions in the medium was determined by 



ICP-MS. The results for Co-ions leaching from the Co- 
Fe NPs showed minimal amount of Co-ions released re- 
gardless to the time incubation or to the presence of 
FCS. 

Effect of Co-Fe NPs on ROS level in Caco-2, TK6, A549, 
MDCK and HepG2 cells 

ROS level in the different cell-lines was explored by the 
DCF method. Table 1 displays an increase in intracellu- 
lar ROS generation in TK6, Caco2 and MDCK cells, but 
not in A549 and HepG2 cells after incubation with 0.5 
and 0.9 mM Co-Fe NPs. The increase is dose-dependent 
and in the high concentration is parallel to the positive 
controls in all three cell-lines; the most sensitive cell- 
line being TK6 which shows 2.4 fold increase of ROS 
after incubation with 0.9 mM Co-Fe NPs. However, in 
A549 cells, which initially seemed very sensitive to 
oxidative stress (ROS increase by 6.3 folds when exposed 
to the positive control of TBHP), their ROS level was 
not affected by Co-Fe NPs at the highest concentration 
(0.9 mM). Similarly there was no ROS increase in 
HepG2 cells. In order to assess the relationship between 
oxidative stress and the averaged toxicity of Co-Fe NPs 
for the five examined cell lines, we determined the linear 
correlation between these two parameters. A linear high 
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Table 1 Fold of ROS increase following incubation with Co-Fe NPs 





TK6 


MDCK 


Caco2 


A549 


HepG2 


Co-Fe NPs (0.9mM) 


*2.43±0.56 


*1.32±0.26 


*1.76±0.15 


# 1.1 ±0.11 


# 0.99±0.16 


Co-Fe NPs (0.5mM) 


*1.86±0.14 


1.21 ±0.21 


# 1.13±0.19 


# 1.01 ±0.07 


# 0.97±0.07 


Positive control TBHP 


2.22±0.48 


1 .46+0.24 


1 .62±0.09 


6.28±0.3 


1.75+0.21 



Caco2, TK6, MDCK, A549 and HepG2 cells were seeded 48 hours before experiments at concentration of 0.2-0.5x1 0 6 cells/well (6 well plate); cells were preloaded 
with 15 uM H 2 DCF-DA for 1.5 h and then exposure for 24 h with 0.5 mM & 0.9 mM Co-Fe NPs; at the end of incubation time cells were washed, trypsinized (if 
adherent), washed and the level of the intracellular fluorescent DCF was analyzed by flow cytometry in terms of of the ratio of the geometric mean of exposed 
cells that of unexposed control cells. The incubation of cells with 0.5 mM TBHP for 45 minutes served as the positive control. The results are the mean of three 
independent experiments ± standard deviation. * Significant differences of each point from the negative control at P<0.05; # Significant differences of each point 
from the positive control at P<0.05. 



correlation (R 2 = 0.97) between the toxicity of Co-Fe 
NPs (0.5 mM, 72h incubation) and the extent of ROS 
generation was observed (Figure 4). A lower linear cor- 
relation (R 2 = 0.68) was obtained at a higher concentra- 
tion (0.9 mM 72 h incubation) of Co-Fe NPs. 

Patterns discovered from the KDD process 

The KDD methodology was applied to the toxicity of Co-Fe 
NPs in the different cellular models. The modeling was ap- 
plied by carrying out training and validation through an it- 
erative process when applying the KDD approach. The first 
KDD goal was to discover rules for determining the toxicity 
of nanoparticles from the experimental results. The input 
data set included the consolidated in-vitro experimental re- 
sults, where each data record has the following attributes: 
(i) Cell type; (ii) Concentration of Co-Fe NPs; (iii) Exposure 
time; (iv) The extent of viability decrease. The model pre- 
dicts the existence or non-existence of toxicity for any com- 
bination of the above parameters. 

The first phase consisted of data cleansing in order to 
validate data correctness, thereby avoiding errors for 
the data-mining process. Data cleansing is part of the 
methodology of data mining. It means that before ap- 
plying the data mining algorithms (such as decision tree 
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Figure 4 Correlation between toxicity and ROS level following 
incubation with Co-Fe NPs. Correlation between the toxicity 
(percent of dead cells) after 72 h incubation with 0.5 mM Co-Fe NPs 
(X axle) versus ROS increase after 1.5 h incubation with 0.5 mM 
Co-Fe NPs (Y axle); the correlation was fitted linearly yielding y = 
0.01 IX + 0.87; all the results represent the ratio of control cells; the 
following symbols represent the different cell-lines: open circle - TK6 
cells, filled circle - Caco2 cells, open triangle - MDCK cells, filled 
triangle - A549 cells, filled rhombus - HepG2 cells. 



here) there are necessary, per project, preprocessing 
stages that look (automatically and by experts) at anom- 
aly and noise in the data, and exclude some data from 
the later data mining step. The step of data cleansing 
causes about 4% of the raw data not to be considered 
due to unreasonable data variation. In the second phase 
we performed data transformation consisting of: (i) 
Normalization of results in terms of control and blank 
results; (ii) Since the results are real numbers while the 
KDD goal is to determine toxicity (binary definition), 
we applied an arbitrary discretization rule where a 
decrease in cell viability, larger than 30% (> EC30) is 
regarded as a toxic response. This level was chosen tak- 
ing into account the measurement error. The validity of 
the threshold chosen was examined by carrying out the 
same modeling with other discretization thresholds of 
25% and 20% for viability decrease. 

The description of the model for toxicity 

The algorithm we applied to model the observed toxicity 
belongs to a type of supervised classifier. We have se- 
lected the J48 decision tree classifier since this classifier 
model can be explained intuitively, in terms of simple if- 
then rules syntax, without any prior knowledge of data 
mining techniques (detailed in the methods). 

In order to evaluate the classifier, we trained in tenfold 
cross-validation mode. This was carried out by splitting 
the data set into ten groups, using nine of the groups for 
training and the tenth for validation, repeating this 
process ten times. This method gives robust result for 
model validation (in essence it does the validation ten 
times, where each time the test set is randomly chosen). 
This known ten-fold cross validation method is consid- 
ered the most robust method for model validation, when 
machine learning is used. This method gives us a better 
estimation of the training error (see Methods), since it 
suffers less from over-fitting. 

The sensitivity and validity of the model 

The sensitivity of the model to the choice of toxicity 
threshold was examined by comparing accuracy and 
kappa coefficient of the model for different thresholds 
(we determine a certain threshold for viability decrease 
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to be defined as toxic or non- toxic). The best fit of the 
model to the data for the threshold values of 30%, 25%, 
and 20% decrease of viability, was performed: The ac- 
curacy of the model, defined as the number of correct 
predictions divided by the total number of measure- 
ments, yielded accuracies of 92.5 % for viability decrease 
threshold larger than 30% (> EC30), 89% for viability 
decrease threshold of 25% (> EC25) and 85.2% for via- 
bility decrease threshold of 20% (> EC20). The kappa 
coefficients are 0.74, 0.67 and 0.61 respectively. Since 
the difference between the three accuracies is signifi- 
cant, the model of viability decrease threshold larger 
than 30% was chosen, since the measurement error 
(SD/mean) was in the range of 5-20%. Due to the low 
level of Co-ion leaching from the NPs, the contribution 
of the Co-ions dissolution was not considered in the ap- 
plied model. 

The validity of the chosen decision tree model was ex- 
amined by comparing it to the well-known "naive bayes" 
classifier. The accuracy of the J48 decision tree model of 
92.5% is significantly higher than the accuracy of 84.2 % 
obtained by the Naive bayes model. The kappa coeffi- 
cient of the Naive Bayes model is 0.48, significantly 
lower than that of the decision tree. 

Table 2 shows the performance matrix obtained from 
the trained decision tree model. Each column of the 
matrix represents the instances in a predicted class 
(toxic or nontoxic), while each row represents the in- 
stances in an actual class. The values of the confusion 
matrix enable computing the classifier performance as 
shown in Table 3. Table 3 also shows the detailed classi- 
fier performance by class, i.e. toxic and non-toxic, 
obtained from the trained decision tree model. The pre- 
cision for a class is the number of TP (true positives, i.e. 
the number of items correctly labeled as belonging to 
the positive class) divided by the total number of ele- 
ments labeled as belonging to the positive class (i.e. the 
sum of true positives and FP (false positives, which are 
items incorrectly labeled as belonging to the class). Re- 
call in this context is defined as the number of true posi- 
tives divided by the total number of elements that actually 
belong to the positive class (i.e. the sum of true positives 
and false negatives, which are items which were not la- 
beled as belonging to the positive class but should have 
been). The F-Measure is the weighted harmonic mean of 
precision and recall {F measure =l/(l/recall+l/precision)). In 



Table 2 Performance matrix 
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b = toxic 
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The confusion matrix obtained from the trained decision tree model. Each 
column of the matrix represents the instances in a predicted class (toxic or 
nontoxic), while each row represents the instances in an actual class. 



Table 3 Analysis erformance measures by class 



TP Rate 


FP Rate 


Precision 


Recall 


F-Measure 


Class 


0.98 


0.3 


0.933 


0.98 


0.955 


non-toxic 


0.7 


0.02 


0.883 


0.7 


0.781 


toxic 



Detailed decision tree classifier performance by class obtained from the 
trained decision tree model. TP (true positives), FP (false positives). 



the current model true positive ratio for the toxic class is 
70% and with F-Measure of above 78%. These are very 
high performance measures, reassuring us of the results 
we have stated in this paper. 

Discussion 

The determination of the extent of Co-ions leaching 
from the NPs due to their dissolution indicates that only 
a minimal amount of Co-ions leached into the medium, 
regardless of the incubation-time or the presence of 
serum proteins. Thus, the highest concentration of 1.2 
mM NPs to which the cells were exposed is expected to 
yield a leaching of 100-110 \iM of Co-ions. Such con- 
centration of Co-ions was shown to be non-toxic con- 
centration towards the various cell models examined in 
this study [14]. Therefore, it is suggested that any toxic 
effect imposed by the NPs should be attributed only to 
the direct effect of Co-Fe NPs themselves. 

The toxicological effects of Co-Fe NPs were examined 
using seven different cell lines, representing various organs 
of the human body. The results show different magnitude 
of decrease in cell viability towards Co-Fe NPs in the dif- 
ferent cell lines as shown by the decision tree model in 
Figure 5. The figure describes the decision tree model 
obtained from the consolidated results after applying the 
KDD process (as described in the results). Applying Deci- 
sion Tree classifier means that the first parameter in the 
tree root (as depicted in Figure 5) is the one, which is the 
most important with respect to the observed toxicity (car- 
ries most of the information about the toxicity). In our 
analysis, the concentration feature becomes the highest 
parameter, as it is the parameter with the maximum infor- 
mation gain, decided by the decision tree algorithm rea- 
soning, with comparison to all other parameters. The 
concentration parameter is divided according to its values. 
Each value is in turn a root of a lower decision tree, recur- 
sively and with the same logic. The model suggests that in- 
cubation with Co-Fe NPs < 200 \xM does not cause any 
toxic effect in all cell types for both incubation times pe- 
riods. Exposure of cells to concentration above 200 \iM 
leads to differential toxicity between the different cell types; 
The toxic effect (more than 30% viability decrease) can be 
seen in lung slices (only 24 h incubation performed), A549, 
NCIH441, MDCK, TK6 and Caco2 cells, however there is 
no toxic effect at all concentrations (200-1200 \iM) on 
dendritic cells and HepG2 cells for both 24 h and 72 h 
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Figure 5 The decision tree model. The decision tree model learned from the consolidated results after applying the KDD process described in 
this work. The decision tree contains different parameters that were used in the data mining procedure: Concentration, Cell-type and time of 
exposure (24 h or 72 h). The outcome/end of each branch is either nontoxic (NT - the light grey squares) or Toxic (T - the dark grey squares). The 
numerical results are given below the outcome (NT or T) in the form of N/n 1; where N represents all data results (NT or T) whereas n ] represents 
the number of data results which do not fulfill the outcome. 



incubation times. TK6 cells seem to be the most sensitive 
cell-line, where Co-Fe-NPs are toxic for TK6 at all concen- 
trations above 200 [iM at both incubation times. NCIH441 
cells are also sensitive to Co-Fe-NPs above 200 \iM but 
only after 72 h incubation. For both MDCK and Caco2 cells 
Co-Fe NPs are toxic at concentration higher than 600 [iM 
after 72 h of incubation. Lung slices as well are sensitive to 
Co-Fe-NPs at concentration higher than 600 [iM after 24 h 
incubation. A549 cells are sensitive to Co-Fe NPs only at 
concentration of 900 |iM and only after 72 h of incubation. 
Thus the cell type turned out to be the second rank influ- 
ential parameter. The third and the lowest rank in the 
model was either the time of exposure (in most cells 24 h 
incubation caused none or minimal effect) or concentration 
depending on the cell type. This model suggests the follow- 
ing hierarchy of cell sensitivity towards viability decrease in- 
duced by Co-Fe NPs: TK6 > Lung slices > NCIH441 > 
Caco-2 = MDCK > A549 > HepG2 = Dendritic. 

The differential viability response of the various tested 
cell-lines towards the Co-Fe NPs may be, in principle, at- 
tributed to their different capacity of interaction with NPs 
in terms of either NP adsorption to the cell surface, NP up- 
take by the cells or both. As the adsorption of Co-Fe NPs 
to the cell surface is non-specific, one would tend to attri- 
bute the observed differential response of the cells to dif- 
ferent active endocytic pathways in the various cell-lines. 
These different pathways can possess dissimilar extent of 
cargo transport as well as different intracellular fate. More- 
over, each cell-line may possess a different repair capacity 



towards the toxic challenge. It should be pointed out that 
even cells originating from the same organ (e.g. NCIH441 
and A549) may retain different viability response. Out of 
the three in-vitro models for lung one would tend to as- 
sume that the lung-slices, which had the highest sensitivity 
among the lung models towards Co-Fe NPs, are the more 
relevant to the in-vivo scenario, when taking into account 
their higher level of resemblance to lung. 

It is intriguing that the two cellular models of the immune 
system (TK6 and primary dendritic cells (DCs)) turned out 
to possess the highest and the lowest sensitivity, respectively. 
DCs are antigen-presenting cells which play a major role in 
the initiation of the specific immune responses, through 
their unique ability to stimulate primary T cells. Immature 
DCs are resident in tissues and epithelial barriers where they 
act as guardians where they detect phagocyte pathogens. 
Such activities imply that DCs are equipped to resist to 
stressing conditions found in inflamed tissues where the 
amount of ROS produced by macrophages or neutrophils, 
for example, is very important. Immature DCs produce fac- 
tors, such as TRANCE (tumor necrosis factor- related 
activation-induced cytokine) and RANK (receptor activa- 
tor of NF-kB), that are associated in rescuing of cells from 
death [16]. Furthermore, DCs express high levels of anti- 
oxidant enzymes, such as Mn-SOD and Prxl, conferring 
the ability to survive in a highly oxidant environment [17]. 
All these features account for their particular survival abil- 
ity to stressing conditions such as incubation in the pres- 
ence of high NPs concentration. The observed much 
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higher sensitivity of lymphocytes vs. dendritic cells is in 
line with previous studies which also have found a high 
sensitivity in TK6 cells; An X-irradiation study found that 
TK6 cells are less efficient in recombinational repair and 
thus have lower resistance to the toxicity of X-irradiation 
relative to WTK1 cells, a human B lymphoblast cell line 
[18]. Another study found that exposure of TK6 to NO, in 
order to imitate the level in inflamed tissue, caused deple- 
tion in reduced glutathione (GSH) level, while the same 
treatment of activated murine macrophages did not inter- 
fere with the GSH homeostasis [19]. 

The hierarchies obtained from the model are an out- 
come of the dose-response curves, where the response 
is an average of viability determined by 2-3 different 
assay methods (MTT, NR and Alamar blue). This hier- 
archy is based on a threshold based on EC30 values. It 
should be stressed that the hierarchy in cell-lines' sensi- 
tivity may depend on the type of NP examined, as well 
as the conditions of incubation; for example gold NPs 
were shown to be toxic to A549 cells but not to HepG2 
and BHK21 (baby hamster kidney) cells [20]. 

Oxidative stress, defined as a situation of an imbalance 
between production of reactive oxygen species (ROS) or 
reactive nitrogen species (RNS) and antioxidant defenses, 
is considered to be an important mechanism of NP- 
induced health effects [21]. It is manifested in the activa- 
tion of ROS, followed by a pro -inflammatory response and 
DNA damage leading to cellular apoptosis and mutagen- 
esis [22] . Therefore, it was important to map the response 
of ROS production of the different cell lines, employed in 
our study. The obtained results suggest that there is a high 
linear correlation (R 2 =0.97) between the toxicity of Co-Fe 
NPs and the extent of ROS generation following their ex- 
posure to Co-Fe NPs (Figure 4). In-vitro studies show that 
NPs generate ROS, deplete endogenous antioxidants, alter 
mitochondrial function and produce oxidative damage to 
lipids and DNA [23]. NPs-induced ROS activation pro- 
motes defense antioxidant response elements. If damage 
proceeds, protective systems are succeeded by mitogen- 
activated protein kinase (MAPK) and NF-i<B-activated 
intracellular signaling, resulting in pro-inflammatory cyto- 
kine, chemokine and matrix metalloproteinase (MMP) re- 
lease leading to apoptosis [22], Our findings provide 
evidence to support the notion that the toxicity of Co-Fe 
NPs is partially due to oxidative stress. 

Conclusions 

In the present study we performed toxicological screening 
of Co-Fe NPs on seven different cell-lines and rat lung 
slices, showing a high correlation to the observed elevation 
of the ROS levels. It was validated that any toxic outcome 
was attributed only to the effect of the Co-Fe NPs them- 
selves (due to a minimal amount of non-toxic Co-ions 



leaching). We were able to integrate and analyze the toxi- 
cological in-vitro data from five different research groups, 
based on averaged threshold values of dose-response 
curves of different viability assays, while maintaining data 
quality control. Data analysis and predictive modeling of 
the obtained data sets were executed by employing a deci- 
sion tree model (J48) where training and validation were 
carried out by an iterative process by applying the Know- 
ledge Discovery from Data approach. The modeling of the 
toxicity data enabled us to obtain a multi-dimensional per- 
spective that cannot be achieved using traditional two or 
even three dimensional plots. It enabled to determine the 
hierarchy pattern of the different parameters studied at a 
threshold of 30% toxicity. Apart from concentration which 
emerged as the expected parameter of the highest rank, 
the findings that the second rank parameter is the cell 
model, as well as the obtained lowest rank for the time of 
exposure were far from being intuitive. Similar findings 
were also recently shown for Co-NPs [14]. The presented 
model agrees with the basic principles of pharmacology 
and toxicology, yet it organizes their relative hierarchy and 
thresholds in a clear model based on the attributes of the 
experimental data. It should be pointed out the present 
study is restricted in terms of its prediction potential to a 
specific type of NPs and to the chosen set of the various 
cell-line. A similar approach, based on decision tree algo- 
rithms, can be applied to different types of NPs. A general 
application of such approach to a much larger database, 
consisting of different NPs and cell-lines is expected to 
provide a more general predictive tool for toxicity assess- 
ment of NPs. In the last case the physical- chemical prop- 
erties of the different NPs will be added as additional 
parameters. Thus, having a very large coherent data-base 
addressing in-vitro toxicity of different NPs employing 
large number of diverse cell lines, taking into account the 
concentration of the particles, duration of exposure and 
the physical-chemical properties of the NPs, we will be in 
a position to provide a general model for predicting the 
toxicity of the NPs. 

Methods 

Materials 

Cobalt chloride, (3- [4,5-dimethylthiazol-2-yl] -2,5- 
diphenyltetrazolium bromide) (MTT), Sodium 
Dodecyl sulphate (SDS), N,N-Dimethyl formamide, 
Neutral Red (NR), Tert-butyl hydroperoxide 
(TBHP), LPS of Escherichia coli 026:B6 were pur- 
chased from Sigma (St. Louis, MO, USA). Dichlorodihy- 
drofluorescein diacetate (Carboxy-H 2 DCF-DA) was 
purchased from Invitrogen (Carlsbad, CA, USA). WST- 
1 was purchased from Roche (Germany). Alamar blue 
assay-kit was purchased from Perkin Elmer. GM-CSF is 
produced by the GM-CSF-transfected J558 cell line (pro- 
vided by D. Gray, University of Edinburg, UK). Murine 
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recombinant IL-6 and FLT-3 ligand (FLT-3 L) were pur- 
chased from Peprotech (Rocky Hill, NJ, USA). The fol- 
lowing antibodies, purchased from BD Pharmingen (San 
Diego, CA, USA) were used: anti-CD40, (3/23), anti- 
TER-119 (TER119) and anti-Grl (RB6-8C5); FITC- 
conjugated anti-I-Ab/d (25.9.17, Ab), PE-conjugated 
anti-CDllc (HL3). Female C57BL/6 mice 3 month old 
were purchased at Charles River (LArbresle, France). 8- 
10 weeks old nulliparous and non-pregnant female 
Wistar Crl:WI (Han) rats were purchased at Charles 
River (Germany). Ferric acetate, (Fe(CH 3 COO) 3 ) was 
purchased from the Shepherd Chemical Company 
(Ohio, USA); Cobalt Acetate (Co(CH 3 COO) 2 .4H 2 0, 
23,7% w/w in Co) was purchased from OMG Kokkola 
Chemicals (Ohio, USA) and diethylene glycol (DEG = O 
(CH 2 CH 2 OH) 2 ) was purchased from Chimica Corona 
(Fiorano Modenese, Italy); All the chemicals were of 
reagent grade and were used without any further 
purification. 

Methods 

Synthesis of Co-Fe NPs 

Co-Fe NPs were synthesized employing the following 
polyol method [24]: cobalt and iron acetates (89,6 and 
179,2 mmol, respectively) were solubilized in 645g of Di- 
ethylene glycol at 110°C for one hour. The solution was 
successively heated to 180°C at a heating rate of 2°C/min 
and then kept at 180°C for three hours. After this grow- 
ing period the dispersion was air cooled to room 
temperature and stored. NPs dispersion is stable over a 
period of two years. Cobalt and iron concentrations of 
the samples were measured after dissolution of the sam- 
ples in concentrated nitric acid, the measurement were 
carried out by solution nebulisation inductively coupled 
plasma emission spectrometry (ICP-AES, Varian Liberty 
51). The final cobalt ferrite concentration was 2.95% (or 
2.86%) (w/v); Where the Co/Fe molar ratio being 0.46; 
29.5 mg/ml cobalt ferrite NPs (2.95% w/v) corresponds 
to a concentration of 4.3xl0 15 NPs/ml with a total area 
of 2.5xl0 12 [im 2 /ml. 

Characterization of Co-Fe NPs 

The crystalline structure of the samples was identified 
from X-Ray diffraction (XRD) patterns recorded in the 
2q range 10-70° with a scan step of 0,05°(2q) for 5s on a 
Philips Xpert pro diffractometer (Cu Ka radiation). The 
crystallite size was determined from the diffraction peak 
by using the Scherrer method. The XRD pattern 
(Figure 6) of the dried sample matches the one expected 
for the spinel phase characteristic of cobalt ferrite. The 
average crystallite size was estimated from the X-ray dif- 
fraction measurement, by using the Scherrer formula; 
we obtain 6.10 nm (average of all the peaks). The 



analysis of the (100) peak gave a diameter of 6.33 nm, 
that could indicate a non-spherical shape. 

DLS measurements were performed on liquids at a con- 
centration of 1 g CoFe204/dm 3 using a Malvern Zetasizer 
nano-S working with a 633 nm laser beam. The measure- 
ment of samples in distilled water (Figure 7) reveals a uni- 
form dispersion of nanoparticles (PDI-0.1) with average 
diameter of 13.6 nm. DLS measurements over a two years 
period did not show any significant alteration of the dis- 
persion stability. 

To investigate size distribution, morphology, agglom- 
eration and the crystal structure in the range of A, we 
use a field emission SEM, the Zeiss 1540 EsB and a con- 
ventional JEOL 3010 operating at 297 kV equipped with a 
LaB6 cathode, post-column Gatan imaging filter and a IK 
slow-scan CCD camera. 

The measurements with the SEM were performed 
using a special STEM detector, which allows studying 
conventional TEM samples and detecting the reflected 
signal above the sample as well as the transmitted sig- 
nal simultaneously. With TEM, the structure informa- 
tion can be obtained by using diffraction patterns that 
allow distinguishing between amorphous and crystal- 
line areas of the specimen and high- resolution images, 
which present the lattice fringes of the crystal. Prelim- 
inary experiments were performed to determine suit- 
able concentrations of the NPs (Figure 8). 

Distribution was also examined by SEM; we chose 
as liquid media DMEM (Dulbeccos Modified Eagle 
Medium) as well as water to measure the influence of 
the physiological medium. All NPs were dispersed in 
deionized water. The state of dispersion was moni- 
tored after some seconds (Figure 9A) or following 72h 
(Figure 9B) at 37°C, 10% C0 2 , vapor saturation or in 
cell culture media (DMEM without supplements), for 
72 h at 37°C (Figure 9C), 10% C0 2 under vapor satur- 
ation, with a dilution of 1:10. We dropped 20 \A of the 
selected NP solution on Lacey-carbon or alternatively 
holey-carbon copper grids laying on filter paper. 

To characterize the crystal structure, transmission elec- 
tron microscopy was used. The zero loss filtered bright 
field TEM image show a high concentration of NPs 
(Figure 10 A). Dark field imaging has yield one particle 
with twinning (Figure 10B). High resolution TEM 
shows mostly single crystals (Figure 10C and D), where 
the determined lattice spacing was 0.26 nm. 

Leaching of cobalt ions from Co-Fe NPs 

The concentration of Co-ions which leached from Co- 
Fe NPs suspensions was measured by Inductively 
Coupled Plasma Mass Spectrometry (ICP-MS, Perkin- 
Elmer SCIEX, Ontario, Canada). We prepared samples 
diluting the stock suspension (120 mM or 28.14 mg/ml) 
in complete culture medium and incubated them for 72 h 
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under standard cell culture conditions (37°C, 5% C0 2 , 
95% humidity) in the presence and absence of 10% FCS. 
After incubation, the separation of Co-ions from NPs 
suspended in culture medium was obtained by filtering 
the samples using centrifugation (3000 rpm corresponding 
to 1810 G for 20 minutes) on Millipore-filter membranes 
of lOkD (Millipore, Italy). In the filtered samples, we mea- 
sured the amount of Co-ions. 

Before the analysis, all the samples were mineral- 
ized by adding the appropriate aliquots of concen- 
trated nitric acid (HN0 3 , ultra pure RS, Carlo Erba 
SpA, Italy), incubated over night at room temperature 
and microwave digested (MDS-2100, CEM Corp. 
USA). 

Cell Cultures and Rat precision-cut slices 
A549 and NCIH441 cells The human alveolar type-II 
(ATII)-like cell lines A549 (ATCC number CCL-185,) 
and NCIH441 (ATCC number HTB-174), purchased 
from LGC-Promochem (Wesel, Germany), were cul- 
tured in RPMI 1640 with L-Glutamine (Invitrogen Cor- 
poration, Germany) supplemented with 10% (v/v) FCS 
(Sigma Aldrich, Germany), 10000 U/ml Penicillin and 
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hydro dinamic diameter (nm) 

Figure 7 Size distribution of Co-Fe NPs. Size distribution of cobalt 
ferrite NPs by DLS. 



10000 U/ml Streptomycin (Invitrogen Corporation, 
Germany). Cells were cultivated in 25 cm 2 tissue culture 
flasks and maintained under standard cell culture condi- 
tions (5% C0 2 , 95% humidity and 37°C in HERAEUS in- 
cubators, Germany). Cells were passaged weekly. 

Caco-2 cells Caco-2 TC7 cells (a kind gift from 
Monique Rousset, INSERM U178, Villejuif, France) 
were cultured in Dulbeccos Modified Eagle Medium 
(DMEM), supplemented with L-glutamine (2 mM), 
20% FCS and 0.05% PSN solution. All cells were grown 
in 75 cm 2 tissue culture flasks (Corning) at 37°C, in a 
humid atmosphere of 5% C0 2 in air. Cells were 
harvested before reaching -60% confluence by using 
0.25% trypsin solution (containing 0.05% EDTA) for 
5 min at 37°C. The cells were centrifuged (1 min at 
400 g, using Sorvall RT6000D), their solution aspirated 
and then re-suspended in growth medium. All culture 
media, antibiotics, trypsin and serum products were 
purchased from Biological Industries (Beit Haemek, 
Israel). 

MDCK and HepG2 cells MDCK and HepG2 cells were 
supplied by the American Type Culture Collection 
(ATCC, USA), mycoplasma-free and with source certifi- 
cation. Experimental cultures were prepared from deep- 
frozen stock vials and always kept in a sub-confluent 
state. They were maintained in complete culture 
medium, prepared for MDCK using DMEM high glu- 
cose (Invitrogen Corporation, Italy) containing 10% v/v 
of Fetal Bovine Australian (Invitrogen Corporation, 
Italy), 2 mM of L-glutamine and 1% v/v of Penicillin/ 
Streptomycin, 10000 U/ml penicillin and 10000 U/ml 
streptomycin (Gibco, Invitrogen Corporation, Italy). For 
HepG2, using DMEM high glucose (Invitrogen Corpor- 
ation, Italy) containing 10% v/v of semi-synthetic FCS II 
(HYCLONE, CELBIO, Milano, Italy) and 1% v/v of Penicil- 
lin/Streptomycin, 10000 U/ml penicillin and 10000 U/ml 
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Figure 8 Characterization of Co-Fe NPs. SEM images of Cobalt- 
Ferrite NPs with different dilutions. 



streptomycin (Gibco, Invitrogen Corporation, Italy). Cell 
preparations were maintained in standard cell culture con- 
ditions (37°C, 5%, C0 2 and 95% humidity, HERAEUS in- 
cubator, Germany) [25]. 

TK6: TK6 cells are human lymphoblastoid cell line 
(purchased from ATCC). The cells were cultured in 
RPMI 1640 supplemented with L-glutamine (2 mM), 
17% FCS and 1% PSN solution. All cells were grown in 
75 cm 2 tissue culture flasks (Corning) at 37°C, in a 
humid atmosphere of 5% C0 2 in air. For each experi- 
ment cells were incubated in 6 wells plates (Griener bio- 
one) at concentration of 0.5 x 10 6 cells/ well. All culture 
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Figure 9 Size distribution of Co-Fe NPs by SEM. Size distribution 
of cobalt ferrite NPs by SEM (A) NPs dispersed in water for 0 h, (B) 
72 h, (C) DMEM for 72 h. 



media, antibiotics, trypsin and serum products were pur- 
chased from Biological Industries (Beit Haemek, Israel). 

Dendritic cells Dendritic cells (DC) were isolated from 
bone marrow cells of C57BL/6 mice as previously de- 
scribed [26,27]. Briefly, bone marrow cells were recov- 
ered by flushing from the femurs. Erythrocytes, TER119 
and GR1 positives cells were removed by magnetic cell 
sorting (Invitrogen, Norway). The remaining negatively 
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Figure 10 Characterization of Co-Fe NPs by TEM TEM images of cobalt ferrite NPs. (A) bright field image, (B) dark field image, (C-D) high 
resolution images, inset- magnification of single particle. 



sorted cells were resuspended at 5xl0 5 cells/ml in 
complete Iscoves modified Dubelccos medium (Gibco, 
Invitrogen, Grand Island, NY, USA) supplemented with 
1% of GM-CSF-transfected J558cell line supernatant, 40 
ng/ml of mouse recombinant FLT3 L and 5 ng/ml of 
mouse recombinant IL-6. On day 3, the cellular super- 
natant was removed and the cells were re-suspended 
under the same conditions. From day 6 to day 11, IL-6 
was removed and FLT3-L was reduced to 20 ng/ml. On 
day 11, the bone marrow cells are differentiated into 
DC, ready for the viability tests. 

Precision cut lung slices Animals (8-10 weeks old nul- 
liparous and non-pregnant female Wistar Crl:WI (Han) 
rats, (Charles River, Germany) were euthanized with an 
overdose of pentobarbital-Na. Lung tissue was prepared 
directly post mortem to conserve the viability of the 



tissue. Through the trachea the lung was carefully filled in 
situ with 10 mL/ 200 g body weight pre- warmed 1.5 % 
agarose - medium solution. The lung was removed and 
put on ice for 20 minutes, allowing the agarose to 
polymerize. Lung lobes were separated, placed on wax, 
and 8 mm tissue cylinders were subsequently prepared. 
The cylinders were placed into the Krumdieck tissue slicer 
(Alabama Research and Development, USA) filled with ice 
cold salt solution (EBSS) and slices with a thickness of ap- 
proximately 200 \im were prepared. The precision cut 
lung slices were then washed three times with pre- 
warmed DMEM/F-12 to completely remove the agarose. 

Cytotoxicity and ROS assays 
MTT assay 

The MTT assay is based on the protocol proposed by 
Dezinot [28]. MTT is a water-soluble tetrazolium salt, 
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which is converted to an insoluble purple formazan by 
cleavage of the tetrazolium ring by succinate dehydro- 
genase within the mitochondria. The cell membrane is 
impermeable to the formazan product and therefore the 
formazan accumulates in healthy cells. The different 
cell-lines were seeded in 96 wells at a non-confluent cell 
density and incubated for 24 h in standard cell culture 
conditions. After 24 h the medium was removed and 
the treatment medium was added with different con- 
centrations of Co-ions or Co-NPs as detailed in the 
results. 

After 24 or 72 h of exposure to Co-Fe NPs, 20 \iL of 
MTT solution (5 mg MTT/lmL PBS) was added to each 
well and the 96 well plates were incubated for 2h at 
37°C. 100 \iL of lysing buffer (containing SDS and N, 
N-Dimethyl formamide pH 4.7) was added to each well, 
and after 1 h cell viability was determined by absorption 
measurement (X=560 nm) and the results expressed as 
the relative viability of cells compared to that of the con- 
trol cells. The Co-Fe NPs' background was subtracted 
from results, where necessary, provided it did not con- 
tribute significantly to the total absorbance/fluorescence 
reading. 

WST- 7 assay for Rat slices 

After an incubation of the slices for 24 h with different 
concentrations of Co-Fe NPs, semi-quantitative determi- 
nations of the mitochondrial metabolic activity was 
performed using the water-soluble tetrazolium WST-1 
assay. The cell viability was determined by absorption 
measurement (\=560 nm) and results were expressed as 
relative viability of cells normalized to that of the 
control. 

Neutral Red assay 

The neutral red (NR) assay employed is based on the 
initial protocol described by Borenfreund and Puerner 
[29], where the accumulation of the neutral red dye in 
the lysosomes of viable, uninjured cells, related to ac- 
tive transport of the dye, is determined. The different 
cell-lines were seeded in 96 well plates at a non- 
confluent cell density and incubated for 24 h under 
standard cell culture conditions. After 24 h the medium 
was removed and the treatment medium was added 
with different concentrations of Co-Fe NPs as detailed 
in the results. After 24 and 72 h of exposure to Co-Fe 
NPs, the treatment solution was replaced with 100 \iL 
of filtered NR stock solution (2 mg NR/1 mL) dissolved 
in appropriate growth medium (1:40) added to each 
well for 2 h in 37°C, followed by rapid washing with 
PBS or medium w/o serum. Then, the cells were agi- 
tated with 100 \iL of a bleaching solution containing 
50% v/v ethanol, 49% v/v ultrapure water and 1% v/v 
acetic acid for 20 minutes. Cell viability was determined 



by absorbance measurement (X=544 nm) and results 
were expressed as relative viability of cells, normalized 
to that of the control cells. The Co-Fe NPs' background 
was subtracted from results, where necessary, provided 
it did not contribute significantly to the total absorb- 
ance/fluorescence reading. 

Alamar Blue assay 

The Alamar Blue assay incorporates a fluoremetric or 
colometric indicator based on detecting metabolic activ- 
ity. Alamar blue (resazurin) is a nontoxic oxidation- 
reduction indicator which is reduced by the metabolic 
activity of the cells reflecting both viable cell number 
and changes in cellular redox activity [30]. The different 
cell-lines were seeded in 96 wells at a non-confluent cell 
density and incubated for 24 h under standard cell cul- 
ture conditions. After 24 h the medium was removed 
and cells were exposed to different concentrations of 
Co-Fe NPs as detailed in the results. After 24 h and 72 h 
of exposure to Co-Fe NPs, 100 \iL of Alamar Blue solu- 
tion (diluted 1:10 from stock solution) was added to 
each well and the 96well plates were incubated for 2-3 h 
in an incubator (37°C). Cell viability was measured in by 
fluorescence (\ e x=485 nm, X em =595 nm) and results 
were expressed as relative viability of cells compared to 
that of control cells. The Co-Fe NPs' background was 
subtracted from results, where necessary, provided it did 
not contribute significantly to the total absorbance/fluor- 
escence reading. 

Determination of intracellular ROS based on 
dichlorodihydrofluorescein diacetate (H 2 DCF-DA) fluorescent 
probe 

Dichlorodihydrofluorescein diacetate (H2DCF-DA) 
(Carboxy-H 2 DCF-DA -Invitrogen D400) was used to 
determine changes in ROS activity. It permeates 
the cell membrane and accumulates mostly in 
the cytosol following deacetylation by esterases 
to dichlorodihydrofluorescein (DCFH). This non- 
fluorescent product is converted by ROS into DCF 
(Ex=488 nm; Em= 525 nm). The different cell-lines 
were seeded in 6 well plates at a non-confluent cell density 
(for the adherent cells) and incubated for 24 h under 
standard cell culture conditions. After 24 h the medium 
was removed and the treatment medium containing 0.5 or 
0.9 mM of Co-Fe NPs was added. After 1.5 h of exposure 
to Co-Fe NPs, the treatment solution was replaced with 
Carboxy-H 2 DCF-DA (15 \xM) by adding it to the free- 
serum medium. After incubation for 1.5 hour, cells were 
washed and the intracellular fluorescence of DCF was 
measured by flow cytometry (FACSCalibur flow cytometer 
from Becton Dickinson Company). The optical interfer- 
ence from the Fe-Co NPs in the flow cytometric measure- 
ments of the dichlorodihydrofluorescein is unlikely as no 
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changes in the side scattering were observed in the pres- 
ence and absence of Co-Fe NPs. Results are expressed in 
terms of Gmean relative to that of control; Exposure to 
TBHP (tert-butyl hydroperoxide) (0.5 mM) for 0.5 hour 
was used as a positive control. 

The statistical significance of the changes for the dif- 
ferent assays of cell viability using specific cell lines was 
determined by two-tailed Students t test. Differences of 
<0.05 were considered statistically significant. 

Methods for knowledge discovery from data 

Knowledge Discovery from Data (KDD) is an automatic, 
exploratory data analysis and modeling of complex data 
sets, such as were available during the present research 
on the toxicological effects of NPs. KDD is an organized 
process of identifying valid, novel, useful, and under- 
standable patterns from large or complex datasets [13]. 
Data mining (DM) is the core of the KDD process, in- 
volving inferring algorithms that explore the data, de- 
velop the model and discover previously unknown 
patterns. The model is used for understanding phenom- 
ena from the data, analysis and prediction. The process 
of KDD starts with determining the KDD goals, and 
ends with the discovered knowledge implementation 
[13]. A brief description of the nine steps of the KDD 
process is detailed in a previous article [14]. 

Data mining algorithms 

Supervised learning methods [13] can be applied to predict 
the toxicity of nanoparticles based on models derived from 
experimental studies. These methods are used to discover a 
functional relationship between independent variables (con- 
centration of Co-Fe NP, cell type, exposure time) and a tar- 
get attribute (toxic or non-toxic), based on a training set, 
which is a set of input (the independent variables) - output 
(the target attribute) pairs which hopefully describe the re- 
lationship in an unbiased manner. 

The decision tree classifier method 

Usually there is a distinction between regression methods 
whose output is real valued (numeric) and classification 
method whose output is categorical (nominal or member 
of an unordered set). In this study we apply the classifica- 
tion methods since we define toxicity in a categorical 
manner (toxic or not- toxic). The specific method chosen 
is based on a decision tree classifier. This is a rooted tree 
(shown in Figure 5) where each of its nodes represents a 
partition of the samples, according to a function of the in- 
put attribute [13]. There are many advantages to using de- 
cision tree for modeling: (1) Decision trees are simple to 
understand and interpret; (2) Decision trees have value 
even with little hard data, important insights can be gener- 
ated; (3) Decision trees can model non-linear relationships 
(4) Decision trees are white box models: if a given result is 



provided by a model, the explanation for the result is eas- 
ily replicated. There are several algorithms for generating 
a decision tree; perhaps the most used one is based on in- 
formation theory concepts: We look for the attribute with 
the highest information gain (concentration or exposure 
time or cell type), and split accordingly, creating a node 
for each of the possibilities, and repeat this process recur- 
sively for each node, with the information gain calculated 
according to the new conditional distribution. There are 
two cases where a node isn't being divided: (1) In all the 
above instances the target attribute is the same; (2) all the 
input attributes are the same, thus not allowing a way to 
distinguish between the different situations; in this case 
we simply create a leaf that predicts the majority. 

The Weka implementation - the J48 decision tree algorithm 

The Weka implementation [31] applied in this work, the 
J48 decision tree algorithm, is based on the well-known 
C4.5 algorithm [32] . This implementation starts with large 
sets of cases belonging to known classes. The cases, de- 
scribed by any mixture of nominal and numeric proper- 
ties, are inspected for patterns that allow the classes to be 
reliably discriminated (i.e. toxic or non-toxic). These pat- 
terns are then expressed as models, in the form of decision 
trees or sets of if-then rules that can be used to classify 
new cases. In order to examine the dependence of the 
obtained results on the model chosen we compared the 
outcome of J48 decision tree algorithm with that of Naive 
Bayes model. The Naive Bayes is a classifier that predicts 
the value of the target attribute, based on probability esti- 
mation [31]. The classifier asks to find the target attribute 
value that maximizes the conditional probability of the 
target attribute given the input attribute. In general, these 
tasks are impractical. However, under the assumption of 
independence of the input attributes (e.g. Co-Fe NPs, con- 
centrations, exposure duration, assay used) for a given tar- 
get attribute (toxicity), these tasks becomes very easy and 
efficient, since the calculation of the desired probability is 
very simple based on Bayes theorem. Finally, as shown in 
the results the J48 decision tree algorithm found to be 
with significantly higher accuracy and higher kappa coeffi- 
cient than the accuracy and kappa coefficient obtained by 
the Naive bayes model, therefore J48 decision tree algo- 
rithm was chosen to be used for our model. We can also 
define the training error (as described below), which is a 
common measure for evaluating the model performance: 
the ratio between the numbers of records wrongly pre- 
dicted after testing the model on the same training set, 
and the whole number of records. 

Empirical estimation of the training error - cross validation 

The main goal of the models described above is to minimize 
the training error [13] - the expected rate of wrong predic- 
tion (unknown since we don't know how the distribution of 
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the population space). Cross-validation is one of several ap- 
proaches to estimating how well the model one has learned 
from some training data [13] is going to perform on future 
as-yet-unseen data. Thus cross validation estimates the 
training error by splitting the known examples into two 
groups: the training set and the test set. First we use the 
training set in order to create a classifier and then we test 
the measure of its misclassification on the test set. This 
method gives us a better estimation of the training error, 
since it suffers less from over-fitting. 
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